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ABSTRACT 

Aims. We present an empirical working model for sunspot umbrae which equally describes observed continuum intensities 
and line profiles. 

Methods. The wings of the infrared Call lines depend sensitively on the temperature gradient at —0.6 < logro.s < +0.3 
but not essentially on the absolute value of T. These lines are observed to remain almost unchanged from photosphere 
to umbra and are thus insensitive to parasitic light. It is also shown that the infrared KI7699A line is suitable for 
umbral spectroscopy since it is not seriously blended, its continuum is well defined and it is less influenced by parasitic 
light as compared to lines in the visible spectrum, due to the smaller umbal contrast. 

Results. Calculations show that the umbral gradient dT/dr, required to fit the Call triplet lines, strongly conflicts with 
the observed profiles of K 1 7699, NaD2 and FeI5434 (g=0), even when assuming vanishing Fell lines for a maximum 
correction of parasitic light. It is shown that the discrepancy from the different line profiles may be removed by adopting 
an opacity enhancement as introduced by Zwaan (1974) from a discussion of continuum contrasts alone. The finally 
proposed umbral working model is very close to a scaled model of the quiet Sun with T e f / = 4000K thus resembling a 
M0 rather than a K5 stellar atmosphere. 
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1. Introduction 

A sensitive method for the determination of the empirical 
model atmospheric of a sunspot umbra is the investigating 
line profiles. The infrared Call triplet lines (at 8498, 8542, 
and 8662A) have so far not been considered. Calculations 
show that the wings of these lines are very sensitive to 
the temperature gradient at —0.6 < logro.s < +0.3 but 
much less sensitive to the absolute value of the tempera- 
ture. Cursory inspection of umbral spectra shows that these 
triplet lines are largely unchanged when compared to the 
spectrum of the neighboring photosphere. This implies an 
almost negligible influence of parasitic light, which anyway 
is reduced with respect to the visible region due to the 
smaller continuum contrast. In addition, the wavelength re- 
gion around 8500A seems not seriously perturbed by purely 
umbral lines. 



2. Observations 

We investigated the infrared Call lines from several 
sunspot spectra observed in summer 1967, 1968, and 1979 
at the Locarno solar station of the Gottingen observatory 
with the plane grating spectrograph (Czerny- Turner type) 
at high dispersion (0.2A/mm) and also with the concave 
grating spectrograph (Pasche-Rungen mounting) at low 
dispersion (2.5A/mm). Microphotometry of the infrared 
Call lines in these umbral spectra consistently shows no 
significant difference with the simultaneously observed line 
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profiles in the neighboring photosphere. For further confir- 
mation of this interesting result we recorded photoelectri- 
cally the line-to-continuum ratio of these Call lines in a 
spot near the disc center on July 7, 1979, and obtain the 
same result. 

Those of our observed Call IR lines which occasion- 
ally show emission cores, always yield symmetric ones 
for Call 8542 but asymmetric ones for 8CaII498 and 
Call 8662 with opposite sign: Call 8498 has a steeper blue, 
Ca II 8662 a steeper red flank of its core emission. 

In addition to these lines, spectra of the strong K 1 7699 
line were investigated in the same manner. In contrast to 
the Call lines, this KI line is not unchanged and there- 
fore strongly affected by parasitic light. A maximum cor- 
rection has been applied by assuming Fell 7711 to vanish 
in umbrae (cf. Stellmacher and Wiehr, 1970). A first pho- 
toelectric recording of this K I line was obtained in the July 
7, 1979, spot with particular care for parasitic light correc- 
tion, which typically consists of 5% penumbral blurring and 
< 1% photospheric scattering (in units of the photospheric 
continuum intensity). 

Figure 1 shows calculated profiles of Call 8542, K 1 7699, 
NaD2 and Fe 1 5434 (g = 0) using different empirical models. 
It can be seen that the strong criterion of the unchanged 
Call lines is only reproduced by Zwaan's (1974) umbral 
model. Closer inspection shows that this agreement is es- 
sentially due to the high temperature gradient in the region 
—0.6 < logTo.% < +0.3, where the line wings are formed. On 
the other hand, it was shown by Stellmacher and Wiehr 
(1975) that Zwaan's (1974) model does not reproduce the 
wings of NaD2 and Fe 5434 within the observational limits 
given by vanishing and, respectively, maximum correction 
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Fig. 1. Observed lines profiles in sunspots: vertical bars give raw and maximum straylight-corrected values. Comparison 
with calculated profiles using the models by Kjeldseth-Moe & Maltby dark (d) and bright (b), Stellmacher & Wiehr (+), 
Henoux (H), Kneer (K), Zwaan with q(X) = 1.05 (O), and model M4 of this paper with q=1.0 (full line). The four panels 
give: [a] Call 8542; [b] K 1 7699 with Zeeman splitting = 3000 G and q=1.3; [c] FeI5434 with q=1.62; [d] NaD 2 (only 
blue wing) with q=1.57. For comparison the photospheric ((dashes) and penumbral ((dash-points) profiles are shown. 



for parasitic light. This also holds when considering our ob- 
servations of the K 7699 line wings (see Fig. 1). The K line 
shows a similar temperature sensitivity as NaD 2 , but it is 
much less influenced by parasitic light due to the almost 
two times smaller umbral continuum contrast at 7700A as 
compared to 5900A. 

We tried to alter the temperature stratification of the 
Stellmacher and Wiehr (1975) model in order to obtain 
a better representation of the infrared Call lines without 
loosing the good fit of this model to the lines Fe 5434, NaD 2 
and also the newly included K7699 line. Our calculations 
show that those temperature stratifications which do fit the 
Ca II lines contradict the observations of either the profiles 
of Fe 5434, NaD 2 and K 7699 or of the relative umbra-to- 
photosphere continuum intensities (of < 0.1 at 5900A and 
< 0.2 at 7700A; see e.g. Zwaan, 1974). Preserving the tem- 
perature gradient required to fit the Ca II lines and shifting 
the T(logr a 5 ) curve towards higher T (at constant logr a 5 ), 
a good representation of Fe5434, NaD 2 , and K7699 is 
achieved. A shift towards lower T, however, is required to fit 
the observed continuum contrasts. Any 'compromise' model 
between these two cases clearly yields Call lines and con- 
tinuum contrasts significantly above, and Fe5434, NaD 2 , 



K 7699 profiles significantly below the range defined by un- 
corrected and maximum corrected observations. 

3. Discussion 

A reasonable fit to all observed lines and continuum data 
(see Fig. 2) could be achieved when introducing the opac- 
ity enhancement q(X) for A < 8500A as suggested by 
Zwaan (1974). Several alterations of the q(X) dependence 
did not improve this fit. The concept of an opacity en- 
hancement for sunspot umbrae had been proposed since 
the wavelength variation of umbral continuum intensities 
at 0.5/i < A < 1.7/1 cannot be reproduced by a model even 
when assuming a maximum temperature gradient from a 
radiative equilibrium stratification. 

On the other hand, Stellmacher and Wiehr (1975) argue 
that the range of observed umbral continuum data and the 
uncertainties of the photospheric reference intensities still 
allow a radiative equilibrium stratification without intro- 
ducing an opacity enhancement. In that case, their model 
M3 gives the 0.5/z intensity at its upper and the 1.7(1 in- 
tensity at its lower observational limits (Stellmacher and 
Wiehr, 1975, cf. Fig. 1). On the other hand, model M3 
nicely describes the observed Fe 5434 and NaD 2 line wings. 
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Fig. 2. Wavelength dependence of the umbra-to- 
photosphere continuum intensity calculated for the 
models M3 (- x - x) and the new model M4 from this paper 
with q(X) = 1.0 (dots), resp., q = q(X) (full line). Bars give 
the range of observations by Ekmann & Maltby (1974, cf. 
Stellmacher & Wiehr, 1975, Fig. 1). 
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Fig. 3. Temperature versus optical depth logr , 5 without 
opacity enhancement for the models by Zwaan (- o - o), 
our former model M3 (- x - x) , and the improved model M4 
from this paper ((full line). 



The consideration of observed infrared Ca II lines in ad- 
dition that of Fe 5434, NaD 2 , K 7699 and the corresponding 
continuum data gives further support to the existence of an 
enhanced opacity, which originally had been proposed by 
Zwaan (1974) for a good representation of umbral contin- 
uum contrasts alone. 

Our present observations do not allow to essentially re- 
fine the q(X) given by Zwaan. In particular, we can not 
verify whether q(X) further increases at A < 5000A, which 
might give some hint to the still unknown physical nature 
of the opacity enhancement. A possible explanation is line 
haze produced by atomic and molecular lines (Zwaan, 1974; 
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Fig. 4. TABLE with characteristic parameters of the final 
umbral model M4; for logro.§ < —0.6 and logro.s > +0.6 
the preent model is identica with our former model M3. 

Gaur et al., 1979) and/or unknown continuum absorbers (as 
e.g. proposed by Oinas, 1974, 1977, for cool stars). 

3.1. Umbral model 

Our improved model 'M4', is given in Fig. 3 together with 
our former model M3 and with Zwaan's (1974) model. The 
new model M4, listed in the table, differs from M3 only 
for the layers —0.6 < logrO.5 < +0.3 in order to fit to the 
observed infrared Call lines discussed in this paper. For 
the outer layers logro.5 < —0.6 M4 is identical with M3 
to preserve the description of observed non-magnetic lines 
as discussed by Stellmacher and Wiehr (1970 and 1972). 
Also the deep layers logro.5 > 0.6 of M4 agree with M3 
in order to meet the radiative equilibrium stratification, 
which is required for logro^ > 0.8 to fit the infrared CI 
lines discussed by Stellmacher and Wiehr (1975). 

The final umbral model M4 then optimally describes 
the whole range of continuum and line profile observations. 
It can thus be considered as an empirical working model 
with wide validity for dark umbral cores. Its temperature 
stratification is close to a scaled model of the undisturbed 
Sun with a scaling factor of T™ b ™ = 4000K and thus to a 
M0 rather than a K5 stellar atmosphere. 
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